
Superfluid qubit systems with ring shaped 
optical lattices 

Luigi Amico\ Davit Aghamalyan^, Filip AuksztoP, Herbert Crepaz^, Rainer Dumke^ & Leong Chuan Kwek^ 

^ CNR-MATIS-IMM & Dipartimento di Fisica e Astronomia, Via S. Sofia 64, 95 1 27 Catania, Italy Centre for Quantum Technologies, 
National University of Singapore, 3 Science Drive 2, Singapore 1 1 7543 and Institute of Advanced Studies, Nanyang Technological 
University, 1 Nanyang Walk, Singapore 63761 6, "^Centre for Quantum Technologies, National University of Singapore, 3 Science 
Drive 2, Singapore 1 1 7543, ^Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 
1 17543 and Division of Physics and Applied Physics, Nanyang Technological University, 21 Nanyang Link, Singapore 637371 , 
^Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 1 1 7543 and National Institute 
of Education and Institute of Advanced Studies, Nanyang Technological University, 1 Nanyang Walk, Singapore 637616. 

We study an experimentally feasible qubit system employing neutral atomic currents. Our system is based 
on bosonic cold atoms trapped in ring-shaped optical lattice potentials. The lattice makes the system strictly 
one dimensional and it provides the infrastructure to realize a tunable ring-ring interaction. Our 
implementation combines the low decoherence rates of neutral cold atoms systems, overcoming single site 
addressing, with the robustness of topologically protected solid state Josephson flux qubits. Characteristic 
fluctuations in the magnetic fields affecting Josephson junction based flux qubits are expected to be 
minimized employing neutral atoms as flux carriers. By breaking the Galilean invariance we demonstrate 
how atomic currents through the lattice provide an implementation of a qubit. This is realized either by 
artificially creating a phase slip in a single ring, or by tunnel coupling of two homogeneous ring lattices. The 
single qubit infrastructure is experimentally investigated with tailored optical potentials. Indeed, we have 
experimentally realized scaled ring-lattice potentials that could host, in principle, ~ 10 of such ring-qubits, 
arranged in a stack configuration, along the laser beam propagation axis. An experimentally viable scheme 
of the two-ring-qubit is discussed, as well. Based on our analysis, we provide protocols to initialize, address, 
and read-out the qubit. 
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A qubit is a two state quantum system that can be coherently manipulated, coupled to its neighbours, and 
measured. Several qubit physical implementations have been proposed in the last decade, all of them 
presenting specific virtues and bottlenecks at different levels^"^. In neutral cold atoms proposals the qubit 
is encoded into well isolated internal atomic states. This allows long coherence times, precise state readout and, in 
principle, scalable quantum registers. However, individual qubit (atom) addressing is a delicate point^ ^ Qubits 
based on Josephson junctions allow fast gate operations and make use of the precision reached by lithography 
techniques^. The decoherence, however, is fast in these systems and it is experimentally challenging to reduce it. 
For charge qubits the main problem arises from dephasing due to background charges in the substrate; flux qubits 
are insensitive to the latter decoherence source, but are influenced by magnetic flux fluctuations due to impaired 
spins proximal to the device^ 

Here we aim at combining the advantages of cold atom and Josephson junction based implementations. The 
basic idea is to use the persistent currents flowing through ring shaped optical lattices^""^^'^^'^^ to realize a cold 
atom analogue of the superconducting flux qubit (see^°'^^"^^ for the different schemes that can be applied to induce 
persistent currents). Recently, superpositions of persistent currents have been thoroughly investigated^^'^^. 

Results 

In this paper we demonstrate how persistent currents flowing in a ring shaped optical lattice can provide a 
physical implementation of a qubit^^'^^'^^. The lattice potential plays an important role in our approach. Indeed, it 
makes strictly one dimensional the atoms' dynamics. Further it provides the means for precise control of the 
confinement and facilitates the qubit-qubit interaction. In our system we break the Galilean invariance. For a 
single ring this is realized by creating a localized 'defect' barrier along a homogeneous lattice^^. Additionally we 
prove that a qubit can be achieved with two homogeneous interacting rings arranged vertically on top of each 
other. In such a system the Galilean invariance is broken along the direction transverse to the two rings. For this 
scheme we analyse the real time dynamics and time-of-flight density distributions. Based on our analysis, we 
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provide viable protocols to initialize, address, and read-out the 
qubits. Indeed, we have experimentally realized scaled ring-lattice 
potentials that could host, in principle, n ~ 10 ring- qubits, arranged 
in a stack configuration, along the laser beam propagation axis. 

Single-ring-qubit: breaking the Galilean invariance on the single 
ring with a site defect. We consider bosonic atoms loaded in a ring- 
shaped potential with identical wells, but with a dimple located at the 
site N — I (see Fig. 1), and pierced by a 'magnetic flux' ^. The system 
is described by the Bose-Hubbard Hamiltonian 

^ i^O i = 0 

where a/s are bosonic operators for atoms trapped in the ring and 
Hi = a]ai. The parameters ti describe the tunnelling between the 
wells along the ring. Since the wells are all identical but one, ti = U 
\ii = 0...N — 2 and - i = t' . Finally, U describes the s-wave 

scattering interactional The 'magnetic flux' is 0= K{z)dz, 

JXi 

where A(z) is the effective vector potential. The effect of the 
dimple is to induce a phase slip at the site N — 1. We assume that 
the density of superfluid is large enough to neglect the fluctuation of 
the number of atoms in each well. In this regime we can assume that 
the system dynamics is characterized by the phases of the superfluid 
order parameter (ff-s, described by the quantum phase modeP^ with 
Josephson coupling// ~ {n)ti {(n) is the average number of bosons in 
each well). The magnetic flux O can be gauged away everywhere but 
at the site {N — l)-th^l Accordingly, the phase difference along 
nearest neighbour sites can be considered small in the 'bulk' and 
the harmonic approximation can be applied. The partition 

function can be written as a path integral: Z= ^ ^[(l)]e~^^^\ where 

the S[(/>] is the Euclidean action. Adapting from the approach 
pursued by Rastelli et al.^^, all the phases except 0 = (f)^_i—<f)Q 



can be integrated out (the integrals are Gaussian). The effective 
action reads 



J 



n2 I 

lU 2{N-1) 



Jo [2U 



(2) 



with the potential U{0) = 



J 



N-1 



{0- -J' cos(^). For large {N - 



I) J' /J and moderate AT, U(0) defines a two-level system. The 
degeneracy point is ^ = 7i: The two states are provided by the 
symmetric and antisymmetric combination of counter- circulating 
currents corresponding to the two minima of U{0). We observe 
that breaking the Galilean invariance of the system provides an 
independent parameter /' facilitating the control of the potential 
landscape. The interaction between 0 and the (harmonic) bulk 
degrees of freedom provides the non local term with G{t) = 
X^^°^^ 7(cL>/)e^^'^, coi being Matsubara frequencies and Y{coi) = 



(N-2)/2 



E 



l + cos[27iA;/(iV-l) 



^ 2/[7(l-cos[27iA;/(iV-l)]) + cof 



The external bath 



vanishes in the thermodynamic limit and the effective action 
reduces to the Caldeira-Leggett one^^. Finally it is worth noting 
that the case of a single junction needs a specific approach but it 
can be demonstrated consistent with Eq.(2). 

Two-rings-qubit: breaking the Galilean invariance with two 
homogeneous coupled rings. We consider bosonic atoms loaded 
in two coupled identical homogeneous rings Fig. 2. We will prove 
that such a system effectively provides a qubit- dynamics (alterna- 
tively to the one-ring qubit implementation discussed above). The 

system is described by the Bose-Hubbard ladder: H = + H^^l + 
Hint, where H^h"^ are the Hamiltonians as in Eq.(l) for the bosons in 
the rings a and b respectively, and 



intensity 
(arb. units) 




100|am 



Figure 1 | Experimental realization of a ring-lattice potential with an 
adjustable weak link (red arrow). Measured intensity distribution with an 
azimuthal lattice spacing of 28 jim and a ring radius of 88 jim (see 
Methods section). The centre peak is the residual zero-order diffraction. 
The effective dynamics of a condensate in such a system is governed by the 
qubit potential as discussed in Eq.(2). The size of the structure is scalable 
and a lower limit is imposed by the diffraction limit of the focusing optics 
(see Methods section). Several rings can be arranged in a stack, along the 
propagation axis of the laser beam (shown in Fig. 5). 
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Figure 2 | Setup for the ring-ring coupling. Two parallel Gaussian laser 
beams (G1,G2) are produced by a combination of two polarizing 
beamsplitter (BSl, BS2). The beam separation D can be controlled by 
moving mirror My. Both beams pass through a lens and interfere to form a 
lattice in z-direction. The distance between the lattice planes is a function 
of l/jy^ which can be varied. The resulting one dimensional lattice is 
combined with vertical beams (LGl, LG2) providing horizontal 
confinement for trapped atoms (See the Methods section). The inset shows 
the ring lattice potentials separated by^^ d = /li^D. The ring-ring 
separation is adjustable by varying the distance D. Such an arrangement 
provides an effective two -level system that can be exploited as a qubit 
(See text). 
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Hint =-g^ (^i^i + • 



(3) 



We observe that along each ring the phase sUps imply twisted 
boundary conditions and therefore they can be localized to a 
specific site, say the N — 1-th. Following a similar procedure as 
employed above, the effective action reads 



a = a,b 

—L—Y^ [rfTrfT'0„(T)G«(T-T')0„(T') 

^ ' n = n h 



(4) 



where each GJ^t) is given by the expression found above for the case 
of a single ring. In this case the phase dynamics is provided by the 
potential 



x = a,b 



2{N-1 



-/cos 



Oa- 



N-2 
N 



with /= {n)g^^. We observe that, for large N, the potential 17(0^, 0^) 
provides that effective phase dynamics of Josephson junctions flux 
qubits realized by Mooji et at. (large A^s corresponds to large 
geometrical inductance of flux qubit devices) In there, the 
landscape was thoroughly analysed. The qubit is made with 
superpositions of the two states and 1^2) corresponding to the 
minima of U{Oa^ Ob). The degeneracy point is achieved by — = 
71. We comment that the ratio /// controls the relative size of the 
energy barriers between minima intra- and minima inter-'unit cells' 
of the (Oa, Oy) phase space, and therefore is important for designing 
the qubit. In our system /// can be fine tuned with the scheme shown 
in Fig. 2. 

Having established that the two tunnel- coupled homogeneous 
rings, indeed, define a two level system, we now study its real-time 
dynamics. We will show that the density of the condensate in the two 
rings can display characteristic oscillations in time. 

We make use of the mean field approximation to analyse the (real 
time) dynamics of the Bose-Hubbard ladder Eqs.(l), (3) (assuming 
that each ring is in a deep superfluid phase). Accordingly Gross- 
Pitaevskii equations are found for the quantities depending on the 
time s. (Pa,i(s) = {(ii(s)) and (pbjis) = (bjis)). Assuming that 

^ ^a,i+i ~^a,i ^^^^ ™g site-independent, we obtain 



— = — Vl—z^ sin 0 

OS 



(6) 



zCOS 0 



where z = (iV^ - NJ/iNa + N^) is the normaUzed 
imbalance between the populations Na and Ni, of the two rings, 
S = 9a — Oi, and s = 2gs{h = l). The parameters are A = 

A^a-M& + ^(^cos^ -cos^^ j g.)^= U/{2g), and p = {Na + 

Nh)/N is the total bosonic density (we included the chemical potential 
^0,). Eqs.(6) can be solved analytically in terms of elliptic func- 
tions^^'^°. Accordingly, the dynamics displays distinct regimes (oscil- 
lating or exponential) as function of the elliptic modulus k, 
depending in turn on A, X, and on the initial population imbalance 
z{0) = Zq. Here we consider the dynamics at Xp^A, i.e. small U/g 
(the analysis of the solutions of the Eqs.(6) in different regimes will be 
presented elsewhere). The results are summarized in Fig. 3. We 



comment that, comparing with A = 0, the oscillations do not average 
to zero (therefore yielding a macroscopic quantum self trapping 
phenomenon^^) and they are faster. The pattern of the circulating 
currents along the two coupled rings can be read out through the 
analysis of the time-of-flight density. As customarily, the spatial 
density distribution in the far field corresponds to the distribution 
in the momentum space at the time when the confinement potential 
is turned off: 



p(k) = 



\w{kx,ky,kz)\^ 



N 



EE E 

/ = 0 j = 0 qe{27in/N} 



kii -xii + 



kii -xii + 



(7) 



N 



2 cos 



k||-X|| +A:,D+( (2+^ 



where wik^MyM^ are Wannier functions (that we considered ident- 
ical for the two rings), k|| -xn = k^ (xi — Xj)-\- ky (ji —yj), Xi = cos(/>p yi 
= sin(/>2 fix the positions of the ring wells in the three dimensional 
space, = 2ni/N being lattice sites along the rings; the expectation 
values involving the Fourier transforms of operators 

aq = lj/VNj2i dindbq = lj/VN are obtained 

for U/t = 0. The density Eq.(7) is displayed in Fig. 4. 

Discussion 

We proposed a construction of flux qubits with atomic neutral cur- 
rents flowing in ring-shaped optical lattice potentials. Persistent cur- 
rents had been experimentally observed in a narrow toroidal trap 
with a weak link^\ The effective action of the system studied in^^ can 
provide a two level system. In contrast with^^"^^, we emphasize how 
we make explicit use of the lattice in our construction, both to confine 
the particles in the rings and to drive the ring-ring interaction. The 
qubits are realized by breaking the Galilean invariance of the system 




time/2g 

Figure 3 | Population imbalance in two coupled rings. We focused on the 
case A^2p. For moderate Zq, oscillations are obtained, with 

co = 2^|a/i + A2 + ^p(zoA-0^) (2A2-l)/[2(l + A2)'/'j| 

corresponding to macroscopic quantum self trapping (blue dashed line). 
The dynamics can be visualized with the h elp of the mechanical system 
provided by a rotator of length ^^^--2^5)^, driven by the external force A. 
The constant solution z(s) = const corresponds to vanishing pendulum 
length (magenta solid line). For A = 0 (inset), the dynamics is 
characterized by Rabi oscillation with Wq = 2g(l -\- kp\/ 1 — Zq / 2^) <co. 
Here kp = 0.1 and A = 4 implying that w ~ Acoq. 
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Figure 4 | Time-of-flight expansion for the two-coupled- rings-qubit. (a,c), vanishing inter-ring tunnelling rate git = 0. In (b,d)> git = 0.9. In the {k„ ky) 
plane the interference fringes with the ring symmetry are due to the momenta of the quantum degenerate gas; the inter-ring tunnelling suppresses the 
interference fringes. In the (ky, k^ plane, ^induces structured interference fringes. The Eq. (7) is calculated for the Bose-Hubbard ladder with 'fluxes' and 
O^, with L7 = 0, and at quantum degeneracy. Results are shown for = 80, = 70, T = O.OS/c^ and AT = 14 with filling fractions of 10 bosons per site. 



either by adding an additional barrier along a single ring lattice 
Eqs.(2), or by tunnel coupling of two homogeneous rings, Eq.(5). 
The latter is proposed to be realized with the scheme in Fig. 2. We 
observe that a suitable variation of such set-up can be exploited also 
to create two qubit gates (each qubit provided by Fig. 1); alternatively, 
a route described in the Methods section can be pursued. 

The analysis of the real time dynamics of such system can be recast 
to a type of coupled Gross-Pitaevskii equations that are characteristic 
for double well potentials, this providing a further proof that the 
system indeed defines a qubit. Accordingly, the basic phenomeno- 
logy of the tunnel-coupled homogeneous rings is demonstrated to be 
characterized by macroscopic quantum self trapping. Since different 
flow states lead to characteristic density patterns in the far field, 
standard expansion of the condensate can be exploited to detect 
the different quantum states of the system (See Fig. 4). 

Our work provides a feasible route to the implementation of a 
functional flux qubit based on persistent atomic currents. For an 
extensive discussion on the one and two qubit gates, please refer to 
the Methods section. The initialization of our qubit can be accomp- 
lished, for example, imparting rotation by exploiting light induced 
torque from Laguerre-Gauss (LG) beams carrying optical angular 
momentum. A two -photon Raman transition between internal 
atomic states can then be used to transfer coherently h. orbital angular 
momentum to the atoms. With this method, transfer efficiencies of 
90% to the rotating state had been demonstrated^ Owing to the 
coherent nature of the Raman process, superpositions of different 
angular momentum states can be prepared^^. Measurements of the 
decay dynamics of a rotating condensate in an optical ring trap 
showed remarkable long lifetimes of the quantized flow states on 



the order of tens of seconds even for high angular momentum (1 = 
10). Phase slips - the dominant decoherence mechanism - condens- 
ate fragmentation and collective excitations which would destroy the 
topologically protected quantum state are strongly suppressed below 
a critical flow velocity. Atom loss in the rotating condensate doesn't 
destroy the state but leads to a slow decrease in the robustness of the 
superfluid where phase slips become more likely^^'^^. 

We comment that, because of the lattice confinement, the gap 
between the two levels of the qubit displays a favourable scaling with 
the number of atoms in the system (assuming that the temperature is 
low enough we can describe the system with Eq.(l))^^'^^'^^. Besides 
making the inter- ring dynamics strictly one dimensional, the lattice 
confinement provides the route to the inter-rings coupling. Indeed, 
the light intensity results to be modulated along the (nearly) cylin- 
drical laser beam. Analysing our experimental configuration, we 
conclude that it is feasible to arrange n ~ 10 ring-qubits in stacks 
configuration (as sketched in Fig. 5) along the beam propagation 
axis. To allow controlled tunnelling between neighbouring lattice 
along the stack, the distance between the ring potentials needs to 
be adjustable in the optical wavelength regime (the schematics in 
Fig. 2 can be employed). A trade-off between high tunnelling rates 
(a necessity for fast gate operations) and an efficient read out and 
addressability of individual stack sites, needs to be analysed. 
Increasing the lattice stack separation after the tunnelling interaction 
has occurred well above the diffraction limit while keeping the atoms 
confined, optical detection and addressing of individual rings 
becomes possible. 

This arrangement produces equal, adjustable ring-ring spacing 
between individual vertical lattice sites and can therefore not readily 
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Figure 5 | Effect of an axial translation on the ring lattice potential, (a) Ring lattice intensity distribution measured at various positions along 
the beam propagation axis around the focal plane (Z = 0). Note that the initial beam, phase modified by the SLM, is not Gaussian any more. The 
optical potential remains undisturbed by a translation of 2.2 times the ring-lattice radius centred around the focal plane (Z = 0). Here R designates the 
ring-lattice radius of 87.5 fim. (b) This is in contrast t o a Gaussian laser beam which exhibits a marked dependence on the axial shift from the focal plane 



where the beam waist co(z) scales with v/ 1 + (2^/2:0) and Rayleigh range Zq. 



be used to couple two two-ring qubits to perform two-qubit 
quantum-gates. The SLM method, however, can be extended to pro- 
duce two ring-lattices in the same horizontal plane, separated by a 
distance larger than the ring diameter. The separation between these 
two adjacent rings can then be programmatically adjusted by updat- 
ing the kinoform to allow tunnelling by mode overlap^^. Combined 
with the adjustable vertical lattice (shown in Fig. 2) this would allow, 
in principle, two-ring qubit stacks to be circumferential tunnel- 
coupled to form two-qubit gates. 

Read out of the angular momentum states can be accomplished 
experimentally with interference of different flow states (i.e. corres- 
ponding to a fragmented superfluid) which maps the phase winding 
into a density modulation that can be measured using time-of-flight 
imaging^^. In the lower panel of Fig. 4 it is shown that different flow 
states lead to characteristic density patterns in the far field. 

We believe that our implementation combines the advantages of 
neutral cold atoms and solid state Josephson junction based flux 
qubits for applications in quantum simulation and computation. 
This promises to exploit the typically low decoherence rates of the 
cold atom systems, overcoming the single site addressing^°, and har- 
ness the full power of macroscopic quantum phenomena in topo- 
logically non trivial systems. The characteristic fluctuations in the 
magnetic fields affecting Josephson junction based fiux qubits are 
expected to be minimized employing neutral atoms as fiux carriers. 

Methods 

Experimental realization of the ring-lattice potential with weak link. We created 
the optical potential with a liquid crystal on silicon spatial light modulator (LC-R 
2500 phase only SLM, Holoeye Photonics AG) which imprints a controlled phase 
onto a collimated laser beam from a 532 nm wavelength diode pumped solid state 
(DPSS) laser. The SLM acts as a programmable phase array and modifies locally the 
phase of an incoming beam. Diffracted light from the computer generated phase 
hologram then forms the desired intensity pattern in the focal plane of an optical 



system (doublet lens, f = 150 mm). The resulting intensity distribution is related to 
the phase distribution of the beam exiting the SLM by Fourier transform. Calculation 
of the required SLM phase pattern (kinoform) has been carried out using an improved 
version of the Mixed-Region-Amplitude-Freedom (MRAF) algorithm^"'^^ with 
angular spectrum propagator. This allows us to simulate numerically the wavefront 
propagation in the optical system without resorting to paraxial approximation. A 
region outside the desired ring lattice pattern (noise region) is dedicated to collect 
unwanted light contributions resulting from the MRAF algorithm's iterative 
optimization process. This can be seen in the measured intensity pattern in Fig. 1 as 
concentric, periodic structures surrounding the ring-lattice and can be filtered out by 
an aperture. 

The ring-lattice potential shown in Fig. 2 and Fig. 5 can be readily scaled down from 
a radius of —90 fim to 5- 10 fim by using a 50 X microscope objective with NA = 0.42 
numerical aperture (Mitutoyo 50 X NIR M-Plan APO) as the focusing optics for the 
SLM beam and with ^2 = 830 nm light, suitable for trapping Rubidium atoms. 
Accounting for the limited reflectivity and diffraction efficiency of the SLM, scat- 
tering into the noise region and losses in the optical system only about 5% of the laser 
light contributes to the optical trapping potential. However this is not a limiting factor 
for small ring-lattice sizes in the tenth of micrometer range as discussed here where 
—50 mW laser power is sufficient to produce well depths of several Free- The gen- 
erated structures are sufficiently smooth, with a measured intensity variation of 4.5% 
rms, to sustain persistent flow-states^\ The barrier height can be dynamically 
modified at a rate up to 50 ms per step, with an upper limit imposed by the frame 
update rate of the SLM LCD panel (60 Hz). 

Setup for the adjustable ring-ring coupling. To allow controlled tunnelling between 
neighbouring lattice stacks the distance between the ring potentials needs to be 
adjustable in the optical wavelength regime. Small distances allow high tunnelling 
rates, a necessity for fast gate operations. This makes it less efficient to read out and 
address individual stack sites, however. Increasing the lattice stack separation after the 
tunnelling interaction has occurred well above the diffraction limit (~/l) while 
keeping the atoms confined, optical detection and addressing of individual rings 
becomes possible. Fig. 2 in the main text illustrates the experimental arrangement to 
produce two adjustable Id ring-lattices by intersecting two Gaussian beams (G1,G2) 
with wavelength /li. The inset in Fig. 2 shows two vertically spaced ring lattice 
potential separated by^*^ d = liflD. The ring-ring separation is controllable by 
changing the beam spacing D between beams Gl and G2, allowing adjustment of the 
ring- ring tunnelling. 

In an experimentally feasible arrangement using light from a Ti:Sa laser at /li ~ 
830 nm, with a beam separation adjustable between D = 10-40 mm and a lens focal 
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length/ = 75 mm, the ring- ring separation can be varied from d = 1.5-6.2 fim. This 
compares to a inter-ring well spacing of 1.5 fim for a ring lattice with 20 lattice sites 
and ring radius of 5 fim. Taking advantage of a large ring-ring separation of 5 fim 
facilitates addressing of individual rings to generate different effective flux-states in a 
stack. Circulation can be created, for instance, with a pulsed pair of Raman beams 
where one of the Raman beams carry H orbital angular momentum. By Raman 
coupling the |F = 2, = 0) and |F = 2,mp = 2) Zeeman ground-states manifolds of 
^^Rh and employing a magnetic gradient field along the vertical axis, the effective two- 
photon Raman detuning can be shifted out of resonance for atoms in rings other than 
the addressed one. The differential Zeeman energy shift between the two Raman 
ground states leads to a magnetic field dependent shift d = fiBgp^fnpB of the two- 
photon Raman detuning. Here /i^ denotes the Bohr magneton, the Lande ^-factor. 
Amp the difference between the magnetic spin-quantum numbers of the two Raman 
states and B the magnetic field strength. With a magnetic field gradient of 180 G/cm - 
a typical value for magnetic traps in BEC experiments - the two-photon Raman 
detuning of a ring which is 5 /im separated from the addressed one with 3 = 0 would 
be shifted by (5 = 126 KHz. As was shown by Wright et aV^, with appropriate choices 
of the magnitude, intensity ratio and detuning of the Raman beams, fractional 
population transfer between the |2,2> <r^ |2,0> states can be accurately controlled by 
varying the two-photon Raman detuning (5 in a range of less than 200 KHz. This was 
demonstrated for Raman beams with Gaussian beam profiles and hence no orbital 
angular momentum was transferred onto the atoms but it can, in principle, be 
adapted for a combination of Gaussian and Laguerre- Gaussian beams to generate 
atomic flux states. 

With a SLM arbitrary optical potentials can be produced in a controlled way only in 
a 2(i-plane - the focal plane of the Fourier transform lens - making it challenging to 
extend and up -scale this scheme to 3d trap arrangements. The experiment, however, 
showed (see Fig. 5) that axially the ring structure potential remains almost undis- 
turbed by a translation along the beam propagation axis of Az = ±2.2 • R, where R 
denotes the ring-lattice radius. The ring-lattice radius is only weakly affected by an 
axial shift along z and scales with ARIR = 0.0097 • z, where z is normalized to the ring- 
lattice radius. For larger axial shifts from the focal plane the quality of the optical 
potential diminishes gradually. Based on our measurements this would allow 
implementation of ring-lattice stacks with more than 10 rings in a vertical arrange- 
ment, assuming a stack separation comparable to the spacing between two adjacent 
lattice sites. Propagation invariant beams may allow a potentially large number of 
rings to be vertically arranged^"^. 

Tunnelling rate estimation for the two coupled ring lattices. The ring lattice 
potential shown in the inset in Fig. 2 can be written as 

Viatt = (/j cosihczf + cosikczf 

Ifpi cos{kiGz)cos{kGz)cos{(j)l)) , 

where/,/ are related to Laguerre functions^". Such a potential with / lattice sites can be 
created directly by diffraction from a SLM or by superposition of two Laguerre- 
Gaussian beams with a positive and negative azimuthal index ± /, respectively^^. The 
WKB estimate of the tunnelling rate gives 



where d = /1//D is the lattice spacing along z-direction. 



(8) 



Demonstration of the one qubit and two qubit unitary gates. The aim of this 
section is to show how the effective phase dynamics of optical ring-lattices with 
impurities serves the construction of one - and two-qubit gates - a necessity for 
universal quantum computation. Here, we adapt results which were obtained by 
Solenov and Mozyrsky*^ for the case of homogeneous rings with impurities. It results 
that a single ring optical lattice with impurity is described by the following effective 
Lagrangian (see Eq. (2) and Supplemental information): 



2U 



■ + 



/ 

N-l 



cos 9 



Then we introduce the canonical momentum P in a usual way: 

dO u 

After performing a Legendre transformation we get the following Hamiltonian: 



(9) 



(10) 



2/1 7' (N-l) 



(0-<I>)^ + cos( 



(11) 



where fi = /7 (7 is an effective mass of the collective particle. The quantization is 

performed by the usual transformation P ^ —dl dO. For 5 = - — — > 1 the 

effective potential in (1 1) can be reduced to a double well; for 0 = 71, the two lowest 
levels of such double well are symmetric and antisymmetric superpositions of the 
states in the left and right wells respectively (See the Supplemental material). The 
effective Hamiltonian can be written as 



(12) 



and the lowest two states are |i/^^> = (0,1)^ and |i/^e> = (1,0)^. An estimate for the gap 
energy can be found employing the WKB approximation^^ 



12^F/U{l-l/5f' 



(13) 



where S > 1. From this formula we can see that the limit of weak barrier and strong 
interactions is most favourable regime to obtain a finite gap between the two energy 
levels of the double level potentiaP^"^"*'^^. We also note that the gap energy splitting can 
be controlled by the height of the impurity barrier. 

Single qubit gates. For the realization of single-qubit rotations, we consider the 
system close to the symmetric double well configuration O ~ ti. In the basis of the two 
level system discussed before the Hamiltonian takes the form: 



O — 71 



(14) 



where {6)qi is the off-diagonal element of the phase- slip in the two-level system basis. 
It is easy to show that spin flip, Hadamard and phase gates can be realized by this 
Hamiltonian. For example, a phase gate can be realized by evolving the state through 
the unitary transformation Uzif^) (tuning the second term of Eq.(14) to zero by 
adjusting the imprinted flux) 



Uz(jS)=exp(/£T(7z) = 



(15) 



After tuning the gap energy close to zero (adjusting the barrier height of the impurity), 
we can realize the following rotation 



LX,(jS)=exp(/aT(7^) = 



cos a z sm a 



I sm a cos a 



(16) 



O — 71 

where a = — - — {0)qiT. When a = 7r/2 and a = 7r/4 the NOT and Hadamard gates are 
respectively realized. 

Two-qubit coupling and gates. The effective dynamics for two coupled qubits, each 
realized as single ring with localized impurity (as in Fig. 1), is governed by the 
Lagrangian 



E- 



!X = a,b 
-j" COS 



2U 



J 



2{N-l) 
N-2 



(0«-O,)'-7'cos(0«) 



N 



(Oa-O,) 



(17) 



Where /' is the Josephson tunnelling energy between two rings. When 0^ = 0/^ = 0 
and j"<^j' the last term reduces to — /" — and the Lagrangian takes the form 



L = J' 



(x = a,b 



1 



V'iN-l) 



]" (ea-e,Y 
~r 2 



(18) 



By applying the same procedure as in the previous section, we obtain the following 
Hamiltonian in the eigen-basis of the two-level systems of rings a and h 



H = Ha + H,+ 'jayMl 



(19) 



(20) 



From this equations it follows that qubit-qubit interactions can be realized using our 
set-up. If we choose the tuning £ ^ 0 and ^^n— -^jj— the natural representation of 
a (SWAPy gate^^ can be obtained: 



U'(t) =exp 



(21) 



where a = — . A CNOT gate can be realized by using two VSWAP gates. It is well 

known that one qubit rotations and a CNOT gate are sufficient to implement a set of 
universal quantum gates'^^ 
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